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Abstract 

0^ ' We consider corrections to the unification of down-quark and charged-lepton Yukawa couplings in 

supcrsymmctric GUTs, which links the large v T — mixing angle to b — > s transitions. These 
^v^j , corrections generically occur in simple grand-unified models with small Higgs representations and 

affect s — > d and b — > d transitions via the mixing of the corresponding right-handed supcrpartners. 
On the basis of a specific SUSY-SO(IO) model, we analyze the constraints from K — K and B<i — B c i 
mixing on the additional dn — sr rotation angle 0. We find that ek already sets a stringent bound on 
6, # max O (1°), indicating a very specific flavor structure of the correction operators. The impact 
of the large neutrino mixings on the unitarity triangle analysis is also briefly discussed, as well as 
(— i | their ability to account for the sizeable CP-violating phase observed recently in B s — > J/tp4> decays. 

1 Introduction 

The start of the LHC at CERN will enable us to study TeV-scale physics directly for the first time. 
CN ■ Most importantly, we will eventually probe the mechanism of electroweak symmetry breaking; moreover, 
we will be able to test various scenarios for new physics beyond the standard model (SM), the leading 
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C > , candidate of which is arguably supersymmetry (SUSY). The presence of supersymmetry at the TeV-scale 
eliminates the quadratically-divergent loop contributions to the Higgs mass and thereby stabilizes the 



o 



electroweak scale against the scales of more fundamental physics. In addition, TeV-scale SUSY models 
provide an attractive mechanism for electroweak symmetry breaking and an appealing candidate for 
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cold dark matter. Furthermore, they offer a compelling outline for the unification of all matter and 
interactions, the first step of which is grand unification Q. 

The near unification of the SM gauge couplings within the minimal supersymmetric standard model 
(MSSM) at M GU T ^ 2 • 10 16 GeV, with the MSSM being valid above the TeV-scale, suggests that the 
standard-model group, Gsm = SU(3) C x SU(2) L x U(l) y , is embedded into a simple gauge group at this 
scale, such as SU(5) §] or SO(10) ||. SO(10) is arguably the most natural GUT group: both the SM 
gauge and matter fields are unified, introducing only one additional matter particle, the right-handed 
neutrino]]] It is an anomaly-free theory and therefore explains the intricate cancellation of the anomalies 
in the standard model [||. Moreover, it contains B — L as a local symmetry, where B and L are baryon 
and lepton number, respectively; the breaking of B — L provides light neutrino masses via the seesaw 
mechanism. Remarkably, Mqtjt is of the right order of magnitude to generate neutrino masses in the 
sub-eV range. Hence, the neutrino masses are linked to the breaking of the GUT symmetry ||. 

Flavor experiments, though not able to access the TeV scale directly, have put strong constraints 
on the MSSM parameters. Due to the lack of deviation with respect to the SM, we expect the new 
sources of flavor mixing and CP violation to be very limited for SUSY particles around the weak scale. 
As formulated by the concept of minimal flavor violation || [7j , we assume that the Yukawa couplings 



1 Strictly speaking, it is the left-handed neutrino singlet. 
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are the only source of flavor violation and (even more) that the supersymmetry breaking parameters 
are universal at some fundamental scale. Within the minimal supergravity (mSUGRA) scenario ||, 
this scale is usually taken to be Mqut- An alternative and arguably more natural choice, however, 
would be the Planck scale, M P1 = G~ N l/2 = 1 • 10 19 GeV.0 The reason to take the MSSM unification 
scale instead is simply that while the use of the renormalization group equations of the MSSM below 
-^GUT is undisputed, the analysis of the region between Mqtjt and Mp\ requires knowledge about the 
grand-unified model. However, the universality of the SUSY-breaking parameters is broken by their 
evolution down to lower energies. Thus the choice of Mgut eliminates potentially important flavor 
effects. In our analysis, we will adopt Mp\ as universality scale, and study consequences of this choice 
in detail. 

In the standard model, fermion mixing is only measurable among the left-handed states and described 
by the quark and lepton mixing matrices, Vckm and Vpmns- Both small and large mixing angles are 
realized: while those in the quark sector are small, two angles in Vpmns turn out to be large. These are 
the neutrino solar and atmospheric mixing angles, where the latter is close to maximal. The effects of 
Vckm and Vpmns are confined to the quark and to the lepton sectors, respectively. In GUTs, however, 
this separation of quark and lepton sector is abrogated as quarks and leptons are unified. Thus their 
masses and mixings are related to each other. While different patterns are possible, it is natural to 
expect imprints of Vpmns on the quark sector as well. In particular, it might be possible to trade off 
small rotations of left-handed down quarks and right-handed leptons against large mixings among right- 
handed down quarks and left-handed leptons, as we will discuss below. The mixing of the right-handed 
fermions is unobservable due to the absence of right-handed flavor-changing currents at the weak scale. 
With weak-scale supersymmetry, the mixing of the corresponding scalar partners of quarks and leptons 
becomes physical. 

The impact of the large atmospheric mixing angle on B s physics has already been investigated in 
detail H [R| [Tl], |l2|. Due to the good agreement of the bottom-quark and tau-lepton masses at Mgut, 
one can adopt the predicted Yukawa unification of down quarks and charged leptons. In order to study 
K and Bd physics, however, one needs to go beyond minimal models and modify the relations among 
the Yukawa couplings]^ Here, we can pursue two avenues: we can either introduce additional Higgs 
fields in larger representations, such as a 45# in SU(5), or parameterize the modifications via higher- 



dimensional operators, suppressed by powers of a more fundamental scale [15, [Tq ]. We opt for the latter 
route for three reasons. One, large Higgs representations introduce a large number of additional fields, 
which both yields large threshold corrections at Mqtjt and makes the gauge coupling blow up shortly 
above the GUT scale. Two, the use of higher-dimensional operators reflects the successful Yukawa 
unification of the third generation; the corrections are suppressed and therefore apply mostly to the 
lighter generations. Finally, we are able to perform a more general study as we do not rely on specific 
Higgs fields. 

In this paper, we will study the impact of the higher-dimensional Yukawa operators on K — K 
and Bd — B^ mixing. A SUSY-SO(IO) GUT with universal super symmetry-breaking parameters at the 
Planck scale will serve as our specific framework. In particular, the precise measurement of ex will 
enable us to tightly constrain the additional (s) quark mixing caused by these operators. The validity 
of our results for more general classes of grand-unified models will also be assessed. 



Alternatively, one might choose the reduced Planck scale, Mpi = (87rGjv) _1,/2 = 2 ■ 10 18 GeV, because it compensates 
for the factor 8ir in the Einstein field equations. 

3 These modifications were neglected in Ref. whose authors consider minimal SU(5). Similarly, Ref. assumes 
a minimal SO(10) model where Vckm describes all SM flavor mixing (the study is from f995, i.e. before the large mixing 
angles in the lepton sector were established). 
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2 Yukawa Unification and Dimension-five Operators 



Grand-unified theories using small Higgs representations to break the electroweak symmetry generically 
predict the unification of down-quark and charged-lepton masses § glf] Before turning to SO(10), let 
us consider minimal SU(5) to bring out the central idea of this work. Here the down-quark singlet, d c , 
and lepton doublet, L, fill up the 5 representation, whereas the quark doublet, Q, as well as the up-quark 
and the electron singlets, u c and e c , are embedded in the 10. As usual, these are lcft-chiral superfields; 
for instance, we have the electron singlet e c L instead of the right-handed electron eR. The adjoint Higgs 
field E breaks SU(5) to the standard-model group, which is then broken to 811(3)^ x U(l) em by a pair 
of quintets, H + H. 

The corresponding Yukawa couplings read 

W Y = Y\ j e ahcde lOf lOf H e + Yi lOf 5 ja H b , (1) 

where a, b, . . . denote SU(5) and i, j flavor indices. The second coupling yields the unification of down- 
quark and charged-lepton Yukawa couplings Yd, e (and thus of the corresponding masses). If Yd, e are 
defined such that the weak doublets are on the left and the singlets on the right, we obtain 

Y d = Yj = Y 2 . (2) 

The mixings of the right-handed (left-handed) down quarks are thus identical (or, more precisely, 
conjugated) to those of the left-handed (right-handed) charged leptons. 

This relation works remarkably well for the third generation but not for the lighter ones. Thus we 
need to include corrections, which are generically generated by higher-dimensional Yukawa operators, 
suppressed by powers of the Planck scale, M P1 Jl|] . With the given particle content, we have two 



operators of mass-dimension five contributing to the down-quark and charged-lepton masses [15] 



Y?i 10f ha ^ H c + ^2 10f he H a . (3) 

The vacuum expectation value (vev) of £ is proportional to hypercharge, (£) = a diag (2, 2, 2; —3, —3) . 
Hence, the second operator modifies the relation (|2|), 

Yd = Y e T + 5^Y CT2 . (4) 

Now we cannot diagonalize both Yukawa matrices simultaneously anymore. In the basis where the 
charged leptons are diagonal, we obtain 

L d D d R d = D e + 5-£- Y ct2 ; (5) 
Mpi 

Dj denote the diagonal Yukawa matrices, and Ld and Rd are unitary rotation matrices for the down- 
quark fields. The good agreement of the bottom and tau masses at the GUT scale indicates that the 
rotation matrices Ld and Rd have a non-trivial 1-2 block only,0 

/* * o\ 

(6) 



Ld, Rd ~ 



* * 

V° V 



The unification of down-quark and charged-lepton masses is a prediction of the SU(4) symmetry, which is present in 
the Pati-Salam model and respected in minimal SU(5). 

f^l ~ Ij it is suppressed with respect to Y| 3 by a /Afpj. 
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Hence, the effect of the additional rotations may only be seen in observables involving the first and 
second generations. 

The effect of the dimension-five operators on proton decay has been studied in great detail [17|. 
In this paper, we point out that the rotation matrix can be severely constrained by the precise 
measurements in K and physics. This, in turn, allows for a complementary study of these operators 
and thus enables us to probe grand-unified models. 

In the following, we will omit the indices of the higher-dimensional operators. For instance, we will 
denote the operators in Eq. (J3j) 

Y? 10, lj -Jj- H = Y£ 10f E ja §^H C + y% lOf E jc ^g- H a . (7) 

Note that these index-less operators represent all possible combinations for the fields to form a singlet, 
and so Y CT is an effective coupling matrix. 



3 Framework 

Let us now turn to SO(10) and consider a model proposed by Chang, Masiero, and Murayama (CMM) 
§. Here the matter fields are unified in the spinor representations, 16^, together with the right-handed 
neutrinos. SO(10) is broken to SU(5) by a pair of Higgs spinors, 16// + 16//. Next, an adjoint field, 
45//, breaks SU(5) and the electroweak symmetry is eventually broken by a pair of fundamental Higgs 
fields, 10// and 10^. In fact, both the SU(5) adjoint and the SU(5) singlet of 45// acquire vevs, the 
latter (denoted by vq ~ 10 17 GeV) being an order of magnitude larger than the former (a ~ 10 16 GeV). 
The Yukawa couplings in the CMM superpotential read 

W Y = 16, Y\ ] 16, 10// + 16, \f 16, + 16, Y^ 16, . (8) 

Let us discuss the individual terms in detail. In the fundamental Higgs field 10//, only the up- type 
Higgs doublet H u acquires a weak-scale vev such that the first term gives masses to the up quarks and 
neutrinos only. The masses for the down quarks and charged leptons are then generated through the 
vev of the down- type Higgs doublet in the second fundamental Higgs field 10'^. (A second Higgs 
field with Yukawa couplings to the SM fermions is generally needed in order to have a non-trivial flavor 
structure.) They are obtained from the second term in Eq. (||) which is of mass-dimension five, in 
contrast to minimal SU(5). As indicated above, this operator actually stands for various, inequivalent 
effective operators with both the SU(5)-singlet and the SU(5)-adjoint vevs of the adjoint Higgs field 45// 
such that the coupling matrix Y2 can only be understood symbolically. The magnitude of this second 
mass term is determined by the vev of the SU(5)-singlet component, vq, which contributes equally 
to down-quark and charged-lepton masses. The strong hierarchy between the t and 6, r masses then 
follows from the vo/Mp\ suppression factor. The smaller SU(5)-breaking vev (a), which is proportional 
to hypercharge as in SU(5), will be important for the modification of the light generation Yukawa 
couplings. The second term in Eq. (@) can be constructed in various ways, for example by integrating 
out SO(10) fields at the Planck scale. The corresponding couplings can be symmetric or antisymmetric 
||18|| ) resulting in an asymmetric effective coupling matrix Y2, as opposed to the symmetric matrices Yi 
and Ytv- Finally, the third term in Eq. (||), again a higher-dimensional operator, generates Majorana 
masses for the right-handed neutrinos. 

We can always choose a basis where one of the Yukawa matrices in Eq. (|8|) is diagonal. In particular, 
the basis where Y^ 7 is diagonal will be referred to as the up-basis. In the CMM model, however, one 
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assumes that and Y^ are simultaneously diagonalizable. This assumption is motivated by the 
observed values for the fermion masses and mixings and might be a result of family symmetries. First, 
we note that the up quarks have a stronger hierarchy than the down quarks, charged leptons, and 
neutrinos. Consequently, the eigenvalues of Y^r must almost have a double hierarchy compared to Yi. 
Then, given the Yukawa couplings in an arbitrary basis, we expect smaller off-diagonal entries in the 
rotation matrices of Yi and Y^v than in Y2 because hierarchical masses generically correspond to small 
mixing. Moreover, the light neutrino mass matrix implies that, barring cancellations, the rotations 
needed to diagonalize Yi should be smaller than those in Vckm [ 19 1 • Thus, even if Y^r is not exactly 



diagonal in the up-basis, the off-diagonal entries in its rotation matrix will be much smaller than the 
entries in Vckm so that they cannot spoil the large mixings among cIr quarks generated by Vpmns- 

Now, with Yi and Y^ being simultaneously diagonal, the flavor structure is (apart from supersym- 
metry-breaking terms, which we will discuss below) fully contained in the remaining coupling, Y2. Let 
us assume for the moment that the relation (|2|) is valid. Then we can rewrite the superpotential in the 
SU(5) basis as 

Wy = 16, D« 16, 10* + 16, (V q *D 2 V e f 16, + 16, D% 16,^f-^ , (9) 

where the second coupling is to be understood as (Q, e c ) T V*D2Ve (d c , L) 45^10^/Mpi (cf. Sec. ||). Then 
V q and Vg coincide with the quark and lepton mixing matrices, Vckm and Vpmns 1 up to phases^] Note 
that the mass matrices of both down quarks and charged leptons have a lopsided structure. 

As discussed in the previous section, the relation (|2|) needs to be modified. Using the SU(5)-breaking 
vev of 45 h, 0") we obtain 

Y d = Yj + 5-y a , (10) 

in accordance with SU(5) discussed above. Again, this notation is symbolic, as Y CT stems from several 
distinct operators. Without these corrections, the large atmospheric mixing angle could directly be 
translated to maximal mixing between the right-handed down squarks 5r and sr. Now the CKM 
matrix diagonalizes Y^Y^ whereas the PMNS matrix diagonalizes Y e Yj, such that we cannot give a 
general relation between the contributions of the correction operators and additional rotations. Let us 
therefore make the ansatz 

R d = (UV e ) T , (11) 

i.e., the rotation of the down-quark singlet fields differs from that of the lepton doublets by a unitary 
matrix U. Clearly, in absence of the correction operators, U = 1. As said before, the goal of this 
paper is to study how much the rotations parameterized by R d differ from those of the charged leptons, 
i.e. whether a sizeable admixture of cIr in s' R is allowed. 

As discussed above, the good bottom-tau unification implies that the (33)-entry of U should be close 
to one, up to a phase, and the remaining entries of the third row and column should be small. Thus we 
parameterize U as 

/cos#e i<?il -sinfle^ 1 -* 2 -^ 3 ) 
= sin6»e^ 2 cosBe^ 3 | , (12) 

y 



U 




6 In the up-basis, Vckm is conventionally denned as the matrix that rotates the left-handed down-quark mass eigenstates 
into the weak eigenbasis, while the inverse of Vpmns rotates the corresponding charged leptons. The transposition between 
Rd and V in Eq. (jH]) is due to relation (^). 



5 



with 9 G [0, 7r/2]. For concreteness, let us assume the tribimaximal form for the leptonic mixing matrix, 
corresponding to the mixing angles 612 = arcsin (l/V3) — 35°, #13 = 0°, and #23 = 45°. In the up-basis, 
we can have V q in its standard parametrization and thereby absorb five of the six phases. Then we can 
indeed identify V q = Vckm- We cannot do so for since we would only move the phases from the 
down-quark Yukawa matrix to the down-squark soft-breaking masses. We therefore choose to have Vt 
with six phases; to see them explicitly, let us write down the mixing matrix for #13 7^ 0, 



/ 



V. 



V 



|ci3 e 4ai 



-7S-73 Sl3e " 
75 "73 Sl3 e " 



73 Cl3e "* 2 
e *(-oi+o 2 +o 4 ) U-±s 13 e 



S13C 



1(6+03) 



\ 



i(-ai+ct2+a 5 ) ( 1 1 



— Cl3 e *( _a i+ Q 3+a4) 

ici3 e i<y ~ ai+a3+aT ^ 



(13) 



where C13 = cos #13 and S13 = sin #13. In this parametrization, we can easily identify the standard phase, 
<5, and then the standard form for Vpmns is given by Vpmns = PlViPr, where 

g i(a 1 -a 2 ) i(cti-a 3 ) " 



Pl 



diag U- ta \ e~ ia \ e~ ta5 



) 



P R = diag 1, e 



(14) 



If ^i3 is indeed zero, the phase 5 drops out of the matrix (|13|). This situation is familiar from the 
standard model: CP violation requires #13 7^ 0. Altogether, for #13 = 0, the mixing matrix for the 
right-handed down quarks in Eq. (|ll|) reads 

/ 2U n e iai - U12 e ia4 2U 2 i e iai - U22 e ia4 e i(-M+o 5 ) \ 



P,i 



1 

71 



V2e ia2 (U u + C/12 e^-^) V2 e ia2 (U 2 i + U22 e i{ - ai ~ ai) ) - V2 e 



73< 



{<t>4 — 01+02+05) 
^4-01+03+05) 



(15) 



with Uij as given in Eq. (12). 

Due to the absence of right-handed multiplets in the standard model, mixing among the right-handed 
down quarks is unobservable. With supersymmetry, however, the mixing of the corresponding squarks 
potentially leads to enhanced amplitudes for flavor-changing processes. As mentioned in the Introduc- 
tion, we will assume universal soft-breaking terms at the Planck scale. This universality, however, is no 
longer present at the electroweak scale. For the scalar masses, this is due to the large Yukawa coupling 
of the third generation in the renormalization group evolution (RGE), such that 

M| (M z ) = diag (m% m% m| (l - A d -)) (16) 

in the case of the cLr soft-breaking terms. The fast RGE between Mpy and vq allows for rather large 
values of Aj 



11]. Now choosing the super-CKM basis where the down quarks are mass eigenstates, 



this matrix is no longer diagonal; in particular, all elements of the 2-3 block are of comparable size: 

/ 1 - sin 2 9 A j/2 sin(20) e _< *« Aj/4 sin 9 e~ i<t>B i Aj/2 \ 

1 - cos 2 9 A.-/2 



m 



\ 



sin(20) e 1 ** A d /4 
A d ~/2 



sm( 



o t( PB r , 



COS ( 



+ Q4 — «5 , 



>} B r , 



sin 

- cos 9 e~ i(t>B = Aj/2 

= <j>l - 4*2 + 4*3 - 4>A + "4 - "5 



(17) 



This observation motivated detailed studies of b — > s transitions in supersymmetric GUT models, in 
particular the decay b — > S7 and B s — B s mixing 10. O]. In the following, we will study the impact 
of the 1-2 and 1-3 blocks, generated by the angle 9 in Eq. (12), on the analogous s — > d and b — > d 
transitions, focussing on K — K and — B^ mixing. 
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4 Meson- Ant imeson Mixing 

The oscillations of a P° — P° meson system can be described by a Schrodinger-type equation, 



dt \\p°(t)) 



i- 



4 ■ (18) 



where M p and T p are two 2x2 hermitian matrices which encode the four transitions P°/P° — ► P°/P° 
via virtual and physical intermediate states, respectively. The physical states and I-P2 ) are obtained 
by diagonalizing M p — The relevant quantity to study new-physics effects in P° — P° mixing is 

the local contribution to the off-diagonal element of M p : 

< 2 =^(m A /= 2 p), (19) 

with Mp, the average meson mass (M Pl +Mp 2 )/2. The effective Hamiltonian H^f- 2 , which comprises 
in general eight effective operators, 

n AF=2 = £ (2Q) 

i=l 

is conveniently expressed at the scale \ip ~ Mp in the and B s systems, and at the scale jjLp < m c in 
the kaon system. For an extensive introduction into the formalism of K — K and s — B^.s mixing, 



see e.g. Ref. [20]. 



One observable which is particularly well-suited to constrain the additional rotation of the dp_ and 
sr squarks in Eq. (|i~2|) is 



which measures the amount of CP-violation in K — K mixing amplitudes. Indeed, \ex\ is very small in 
the standard model and its experimental value, measured with high precision, leaves only little room 
for new physics. The correction factor K e above parameterizes both the small deviation of sm.<j) eK = 
AM K /(AMj, + Ar^/4) 1 / 2 from l/S/2 and the small contribution from the phase of the isospin-zero 
K — ► 7T7T decay amplitude. This factor was estimated to k 6 = 0.92 ± 0.02 [SlJ assuming the standard 
model. Its modification in the presence of new physics will not alter our analysis, and we will ignore 
this complication. The mass difference AMk between the two eigenstates Kl and K$ receives both 
short-distance and long-distance contributions, such that the constraint on possible new-physics effects 
in the short-distance part, 

(AM K ) SD = 2Re{Mf 2 ), (22) 

is somewhat diluted among hadronic uncertainties. Despite its precise experimental knowledge, AMk 
will thus play a minor role in our study. 

On the contrary, when new sources of CP-violation in the kaon system are small, two observables 
in the B^ system will prove useful to gain information on the mixing angle 6. These are the mass 
difference, 

AM d = 2 M Pd , (23) 
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AM^ xp = (3.483 ± 0.006) • 10~ 12 MeV §| 
AMf p = (3.337 ± 0.033) • 10" 10 MeV || 
AM s exp = (117.0 ± 0.8) • 10- 10 MeV 



| e ^|ex P = (2.229 ±0.012) • 10' 

S e *P =0.671 ±0.024 

j/i)K s 

,exp _ ( n ^^ + 0.29 



in 



(-0.77lg;g)U(-2.36lg;g)rad 



Table 1: Current experimental values of the various AF = 2 observables considered in Eqs. (21- 2< 



and the coefficient of the sin (AM^t) term in the B& — > J/ipKs time-dependent CP asymmetry, 



S 



JI^K S 



sin (2/3 + <f>£) ~ Im 



1 M$ 1 


, /3 = arg 


f ^1 


V 


Mf 2 d 


y 


l v* cd v ch \ 



(24) 



The phase 4> d parameterizes CP-violating effects beyond the SM in B d — B d mixing. Here and in the 
following, we use the standard CKM phase convention. 

Finally, we will also consider the mass difference in the B s system, 



AM, 



M 



(25) 



as well as the phase measured in the B s —* J ftp cj) time-dependent angular distribution, 



-2@f = -2/3, + <f>f ~ arcsin Im 



M 



12 



M 



12 



0s 



arg 



v*v cb 



1VI 12 



In the SM, f3 s is tiny: 2/3 s ~ 0.04. As long as is not too small, we thus have — 2/3| ff ~ 4>^. On 
the other hand, one also has <p s = arg(— M^ s /rf 2 s ) ~ <j)f [p2fl . In the following, we will thus identify 

<j> a = -2/3f . 

The current experimental values of the various observables above are reported in Tab. [l]. 
4.1 Standard-Model Contributions 

In the standard model, W box diagrams with virtual t and/or c flavors generate the effective operators 
Qk LL = (d Ll „s L ) (d L -fa L ) , Q^L = { q L%ibL) (q Ll » bL ) (27) 



for kaons (see Fig. |l(a)| ) and B q (q = s or d), respectively. The corresponding Wilson coefficients at the 
scale [ip read 

CX LL (hk) = Wk^k) [(V c * d V cs ) 2 m S (x c ) + 2(V c * d V cs )(V t * d V ts ) m S (x c ,xt) + (V t * d V ts ) 2 m So(x t )] , 
C^ L (VB q ) = AUB^B^V^fnBS^xt), 

where the factors 

-2/9 



U K {ji)= afXli) 



and = CifXlA 



-6/23 



u A,' 



(29) 



encode the fj,jc, fj>B q -dependent parts of the short-distance QCD corrections up to next-to- leading order 
(NLO), while r/j account for their fj,x, fJ-B q -independent contributions [25, 26]; their values are given in 
Tab. |2[ The loop functions So(x q ) and So(x c ,xt) are listed in the appendix. Finally, x q = m q /My? and 
m q = fn q {m q ) is the MS mass. 
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(a) 

Figure 1: Dominant short-distance contributions to Mf^ (a) in the SM; (b) in the CMM extension. 



In order to compute M^ Bq , we still need the matrix elements of Q^ L and Qp]^ L - These are 
parameterized in terms of "bag factors" Bp, computed at the scale fi = O(fip): 



p°\Q Y p Lh (ri\P°) = \m*fIb p {u.i 



(30) 



where Fp is the decay constant of the P meson. The scale dependences of Up(fi) and Bp(fi) cancel 
each other, so that it is convenient to define the renormalization-group- invariant parameters Bp = 
Bp(n)U P (jj,). Eqs. d|), d|), and @) then lead to 



(31) 



n 2 m 2 

{m k )S m = _E_3L MkF 2 kBk [(aSJVSoCsc) + 2(X c ds )(\ t ds ) m So(x c ,x t ) + (X ds ) 2 V2 S (x t )] , 



(< 2 9 ) SM = ^^M Bq Fl q B Bq {X\ b fn B S,{x t l 
where one defines A^. = V^Vkj- 

4.2 CMM Contributions 

In the context of the CMM model, the dominant supersymmetric effects originate from gluino box 



diagrams with virtual dp, sr, and bp flavors due to the large mixings in Eq. (15). This gives rise to the 
parity-reflected operators (Fig. 1(b)) 



Qk RR = (dR7»s R ) (dnrfsn) , Qp[ RR = (q R ^b R ) (q R ^b R ) . (32) 
The initial conditions for the Wilson coefficients at the SUSY scale Ms = 0(m d .,mg) read |9|, [PH 



L K,B a \ M S) 



16vr 2 a 2 s (M s 



^ ( R d)mj(Rd)mi(Rd)kj(Rd)ti L o( r m,r k ) 



(33) 



where = (1,2) in the kaon case, (1,3) in the B d case, and (2,3) in the B s case. The loop function 
Lo(r m ,rk) is defined in the appendix, the down- type squark mixing matrix R d was given in Eq. (H), 



and rj = m ~ /m~. Exploiting the mass degeneracy of the first two generations (see Eq. fll6|) ) as well as 



the unitarity of R d , Eq. fl33|) simplifies to 
Wtt 2 a 2 s (M s 



C K,B a \ M S) 



G 2 F M& 2m? 



[(R d ) 3j (Rd)h] {Loin, n) - 2L (n, r 3 ) + L (r 3 , r 3 )} , 



n = m\j m 2 g, r 3 = m|(l - A j) / m~ . 



(34) 
(35) 



The RGE of the above Wilson coefficients from the scale Ms down to the scale fiK,B q is performed in 
two steps: first, the leading-order matching coefficients in Eq. ( |34|) are evolved down to fit = 0(mt) by 
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means of the leading-order RGE factor rjQ = [a>s (Ms) / ctg (fJ't)] 2 ^ 7 ■ The remaining evolution, running 
over two orders of magnitude, is achieved using NLO formulas - essentially the Uk{hk), V2, UB q {^B q )i 
and t]b factors of Sec. 4.1. The 0(a s ) QCD corrections to the SM function So(xt) at the scale fit, which 
are contained in 772 and rjB, should be removed. Denoting them by r = 0.985 |26[j, we get 



C% MM M = UkM m - % C£ MM (M 5 ), (36) 

r 

and similarly for C^^ M (^s 9 ). The cancellation of the /^-dependence between the two parts of the 
evolution is of course incomplete, yet this is a numerically small effect which can be neglected. 

The bag parameters of the effective operators Qx RR and Q)} RR are identical to those of the SM 
operators in Eq. (|30| ) such that the CMM contributions to the matrix elements Mf 2 finally read 

-2i(f>f 



(MS™ = ^M K FlB K 6 .f' 2 *' 2S* S«) (ri , r3 ), 
6m~ lb r 

{M s. r = ^Sl MBA B Bd £!!f^ 35* S « (ri , r3) , (37) 



6m~ 



(Mif) oMM = ^ MbM b b , e- 2 '*-™ 1 " asm s ® (n , r3) , 



6m? 



where we explicitly display the factors (Rd)si in Eq. (|4|), S^'(ri,rs) = Lo(n,ri) — 2Lo(ri,r3) + 
Lo(f3,r^), and the CMM phases (f>K, 4>B d , and <pB s have been defined in Eq. (|l7|). Note that they fulfill 
the relation 4>b h = 4>K + <Ab s - 



4.3 Additional Supersymmetric Contributions 

Finally, we comment on the supersymmetric contributions which do not exhibit the large enhancement 
factors characteristic of the CMM model, namely charged-Higgs(-ff)-quark and chargino(x)-squark box 
diagrams. They do not introduce new operators, and the flavor structure of the corresponding matrix 
elements is the same as in the SM, 

C 2 M 2 

(M«) H+ * = -^-fM K F 2 K B K {2(X c M s ) V i S H (c,t) + (X ds ) 2 m [S H (t,t) + S x (t,t)]} , 

r 2 m 2 ^ 

(M%) H+ * = ^^M Bq FlB Bq (Xl b ) 2 r) B [S H (t,t) + S x (t,t)] . 

The loop functions Sh(c, t), t), and S x (t, t) are given explicitly in Ref. §. The factor jyf = 0.21 § 

denotes leading-order QCD corrections to the charged-Higgs box with virtual flavors (c, t). Numerically, 
charged-Higgs and chargino contributions are small compared to CMM effects. We checked explicitly 
that they can be neglected in our analysis. 



5 Numerical Analysis 

We are now ready to investigate the constraints of K — K and — B^ mixing on the angle 6 in the 
down-type squark mixing matrix R^. Since we do not expect a miraculous cancellation of the phases 
4>i and 02, we will first focus on the case where sin2<^ ~ 0(1) (Sec. |5.2| ) and derive constraints on 
from \ex\ alone. We will then turn to the special case sin2^^- ~ (Sec. |5.3| ) where, as we will see, 
interesting constraints can still be obtained from AMk, AMj, Sj/^k s ^ and AM^/AM S . 
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K e = 0.92 ± 0.02 

F K = (156.1 ±0.8) MeV 

B K = 0.75 ± 0.07 

(270 ± 30) MeV 



— 



(1.21 ±0.04) 



m c {m c ) = (1.266 ±0.014) GeV 
fn t (m t ) = (162.1 ± 1.2) GeV 
a s (M z ) = 0.1176 ±0.0020 



II 



\V US \ = 0.2246 ± 0.0012 
jV^I = (41.6 ±0.6) • 10" 3 
\V ub \ = (3.95 ±0.35) • 10~ 3 
7 = (70.7; 



. + 5.7\o 
■ 7.0) 



•qi = (1.32 ±0.32) 

T] 2 = 0.57 ±0.01 
r] 3 = 0.47 ± 0.05 



L30 GeV 

m c (m c ) 



1.1 



VB 



0.551 ±0.007 



3| 

see text] 



Table 2: Input parameters. 



The values of the various input parameters adopted in our numerical analysis are reported in Tab. 
Inputs related to CKM elements have to be protected from new-physics impact. To this end, we 
determine the CKM matrix from the elements \V ub \, \V c b\, \V US \, and 5, the CP-phase in the standard 
parametrization, which equals the angle 7 of the unitarity triangle to very good accuracy. The three 
CKM elements are extracted from tree-level decays. We use \V US \ = 0.2246 ± 0.0012 [27], the inclusive 
determination \ V cb \ = (41.6 ± 0.6) • 10~ 3 f2|, and the average of inclusive and exclusive determinations 
\V ub \ = (3.95 ± 0.35) • 10~ 3 || . The angle 7 is determined via 7 = vr - a - /? = vr - a eff - /3 cfr , 
with p eS = (3 + <j>%/2 = (21.1 ± 0.9)° from Sj /i>Ks @ and a cS = a - <j>f/2 = (88.2±2;|)° from 
B — > 7T7T, 7T/3, pp decays p3]. The dependence on the new-physics phase <fi^ cancels out in the sum 
a cS + (3 cS , such that 7 = (70.71 7'J) is indeed free from new-physics contamination. 

No assumption is made on the squark mixing parameters 9, (fix, ^B d ^ and <j)B a prior to the analysis of 
the observables in Tab. [I]. The supersymmetric parameters (in particular nig, n, and or equivalently 
m g , rrij, and Aj), on the other hand, are chosen such as to satisfy the constraints coming from other 
observables. The identification of viable sets of SUSY parameters is the subject of the next section. 



5.1 CMM Parameter Sets 

In the CMM model, the large number of free SUSY parameters shrinks to six input parameters at the 
electroweak scale (in addition to 6 and the CMM phases (fix, 4>B d , and 4>B a )- These can be chosen as the 
gluino mass m g , the first-generation d,R and ur soft masses mj and the ratio of the (ll)-elements 
of the trilinear and Yukawa couplings in the super-CKM basis a\ = (Adjii/ '(Yd)n, the phase of the 
p parameter in the Higgs potential arg(/x), and the ratio of the two Higgs-doublet vevs tan/3. The 
RGE links these CMM inputs to the remaining SUSY parameters via the assumption of universal soft- 
breaking parameters at the Planck scale and the intermediate SO(10) and SU(5) GUT relations. Note 
that the similar input parameters in the CMM model and in specific SUSY scenarios without grand 
unification still lead to very different phenomenologies. In such well-studied scenarios as mSUGRA or 
the CMSSM, the SUSY-breaking parameters are universal at Mgut, as mentioned in the Introduction, 
leaving the universal gaugino and scalar masses, and mo, the trilinear coupling A, as well as the 
sign of p and tan j3 as free parameters. In contrast to GUT models, however, these scenarios do not 
relate quarks and leptons to each other; the MSSM fields can be rotated independently and the large 
lepton mixing angles do not become visible in the quark sector. 

To establish benchmarks for our analysis of the down-squark mixing angle 9 in K — K and — B^ 
mixing, we make sure that the chosen CMM input parameters are in accord with the other observables 

7 The specification of both mj and ma fixes the D-term scalar mass splitting [(L2L |35| . 
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1000 ■ 

300 400 500 600 700 800 900 1000 

Figure 2: Down-squark mass splitting Aj as a function of m g and [GeV]. White: negative soft 
masses. Black: excluded by lower bound on light Higgs mass. 



sensitive to CMM effects, and that they respect constraints common to generic SUSY scenarios. To this 
end, we make use of the Mathematica code written by the authors of Ref. [12|, which implements the 
relations between the CMM input parameters discussed above and the remaining SUSY parameters at 
the electroweak scale. The most restrictive observable is the experimental lower bound on the mass of 
the lightest Higgs boson m^. For small values of tan/3 it is close to the SM bound, m-h > 114.4 GeV 
||36| . The main radiative corrections to the tree-level Higgs mass in the MSSM, m^ rcc < Mz|cos2/3|, 
stem from (s)top loops. For very small values of tan/? ~ 3 the large top Yukawa coupling in the RGE 
drives the stop mass to low values, such that the Higgs mass bound cannot be fulfilled. In our analysis 
we choose tan /3 = 5, such that the top Yukawa coupling gets smaller, but the natural hierarchy between 
the top and bottom Yukawa couplings, induced by vq/Mp\ in the CMM superpotential, is preserved. We 
fix the inputs orjmg = 1.8 and arg(/i) = 0, such that the allowed space for m g and mg around 1 TeV 
is large. Finally, we take = mg as in Ref. [12]. In Fig. || we show the mass splitting parameter Ag 
in the m g — mg plane for this scenario. Black regions are excluded by the Higgs mass bound. White 
regions are forbidden due to negative soft mass parameters. Additional constraints arise from processes 
reflecting the large atmospheric neutrino mixing angle like r — > fij, b — > 57, and the mass difference 
AM S ; these can cut further into the low m g and mg regions. 

Based on these considerations, we select three sets of CMM input parameters, given in Tab. |. As 
said above, these parameters are defined at the electroweak scale, more precisely at Mz, in Ref. [12|. 
For consistency, we will thus set M5 = Mz (and correspondingly r?6 = 1, neglecting the small effect 
of mt 7^ Mz) in our analysis of meson- antimeson mixing. Sets 2 and 3 do satisfy the AM, constraint 
for all values of 6 and while Set 1 requires |2<^b s | to be between 1.2 and 2.4 radians for small 8 
to satisfy this constraint. Note that especially Set 1 (with small m g and large A;) is chosen such that 
CMM effects in b — > s, b — > d, and s — > d transitions are large. 
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m-g [GeV] 


m d ~ [GeV] 




jmax To 1 


Set 1 


400 


2000 


0.52 


0.5 


Set 2 


700 


2000 


0.44 


0.9 


Set 3 


700 


3000 


0.51 


0.9 



Table 3: CMM parameter sets for fixed a d /m^ = 1.8, arg(^) = 0, and tan/3 = 5, satisfying the 
constraints discussed in Sec. 5J_. The last column shows the maximal mixing angle Q max allowed by \ex\ 
for sin 2^ = 1 (the symmetric solution 9 £ [(it — # max )/2, it/ 2] is excluded by B physics observables, 
see Fig. ||). 



5.2 Scenario I: sin 2(f) K ~ 0(1) 

As long as the CMM phase cftx is n °t too close to zero, \ex\ gives the best constraint on 9. The 
dependence of # max on the relevant combinations of parameters, i.e., sin 2^ /m| , m^/m g , and Aj, is 
summarized in Fig. ^-left. The plain black and dashed gray lines (which happen to be nearly superposed) 
correspond to mg/m g and Aj of Set 2 and Set 1, respectively, while the two other lines are obtained by 
interchanging m^/m g . As one can see, for | sin 20^-|/m? > 1 TeV -2 and typical values of the parameters 
m dl m 9 an< ^ # max is of the order of one degree. Fig. || has been obtained treating the errors in Tab. [2] 
as flat, yet a different error treatment - and/or inflated errors in Tab. |2|- would not change this picture 
significantly. Fixing (fix to vr/4, the precise limits obtained for the various parameter sets defined in 
Sec. [54] are displayed in the last column of Tab. ||. The small contributions in Sec. 43 have no impact 
on these numbers. 

In the B d and B s systems, the SM contributions are not as suppressed as for ek- Consequently, the 
smallness of (9 max prevents any visible effect in AM d and Sj/^,x s , while the formulas for AM S and 4> s 
are well approximated setting 9 = 0. Interestingly, sizeable CMM contributions in the B s system may 
be welcome to reduce the 2.2u discrepancy between the SM prediction for <p s and its experimental value 
[24]. Within Set 1 it is possible to bring this discrepancy down to the one-sigma level while satisfying 
all existing constraints, see Fig. fright. 

Finally, we briefly comment on the dependence of # max on the hypothesis of tribimaximal lepton 
mixing. In particular, one might expect the 23-mixing angle to be large but not it/ 4. In this case, 
Im [(Rd) 32 (Rd) 3 i\ 2 = — I sin 4 023 sin 2 (29) sin(20#-) for #13 = 0. Hence, for large #23, the constraints on 
9 do not differ much. For a sizeable 13-mixing angle in Vi, \ex\ gets additional contributions: 



A (im [(R d ) 32 (R d )* 31 ] 2 ^j = sm9 13 sin 3 9 23 sin (29) [sin(2^) cos(2fl) cos(0 3 - 2 + «4 - <*i - S) 



cos(24>k) sin( 



+ a 4 - ai - 6)] + O (sin 2 9 13 ) . 



(39) 



No large numerical factors offset the sin #i3-suppression, so that the modified 9 bounds are again as 
stringent as those exemplified in Fig. ||[ 

Up to now, we have taken the viewpoint of a fixed sparticle spectrum, and investigated the correlation 
between effects in b — > s and b, s — > d transitions governed by the mixing angle 9. As 9 turns out to be 
restricted to very small values, it is interesting to consider the opposite viewpoint of a fixed 'natural' 9 
value - say, sin# = 0.5 - and derive the corresponding constraints on sparticle masses from ex- Setting 
again (f>x to vr/4, we find that a soft mass scale m g ~ 2TeV is possible only if the ratio m^/m g ~ 1. 
In such a scenario, however, the mass splitting parameter Aj is very small (cf. Fig. ||), such that 
CMM effects in other observables are negligible. For larger values of the ratio m^/m g , Aj increases and 
accordingly the constraints on m g are much more stringent (for example m g > 20TeV for m^/m g = 2). 
CMM effects in B d and B s physics are thus again killed, this time by the strong l/m~ suppression factor. 
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sin20 K /m? [TeV" 2 ] 20 Bs 

Figure 3: Left: Constraints on 9 from \eic\ as a function of sin 2<j> K /m| for A d ~ = 0.44 (black) / 0.52 
(gray) and m^/wig = 2.86 (plain) / 5 (dashed). Right: CP-violating phase 4> s as a function of the 
CMM phase 4>b s - The light gray (dark gray) curve corresponds to Set 1 (Set 2) with 9 = 0. The SM 
prediction (horizontal line) is recovered for 2(f>B 3 = 0, ±tt. The broad gray band indicates the one-sigma 
measurement / [2^/ . 



5.3 Scenario II: sin2</>K: ~ 

If sin20i^ is close to zero, CMM effects cannot make their way into ImM^ anymore, and the best 
constraints on 9 are obtained from AMk and B physics observables. As mentioned in Sec. ||, AMk 
is plagued by hadronic uncertainties, so that we merely impose |AM^ MM | < AM^ xp to stay on the 
conservative side. In this case, for m g ~ 700 GeV, the constraint from AMk only starts to compete 
with that from \ek\ when \(Pk\ = O(0.1°), corresponding to # max ~ 10° — 30° (depending on the precise 
values of Aj and m^/wig). The constraints from AM^, Sj/^k s , and AM^/ AM S are in general better, 
as we illustrate in Fig. |I| for Set 1 and Set 2. Note that the constraint from AM^/ AM S depends on 
both (f>B s and 4>B d = <\>K + §B S - The plots shown in Fig. [| correspond to 4>k = and 4>k = k/2. Other 
<pK values lead to different plots, with however the same general appearance, in particular the exclusion 
of small 9 angles for some specific <pB s values. For these specific values, the tight bounds on 9 derived 
in Sec. are thus even surpassed. 

As mentioned previously, <j) s can cut further into the parameter space, especially for negative (pB s 
values, see Fig. ^-right. However, this does not change the typical value of # max obtained from B physics 
observables, which is of ten or a few tens of degrees. 



5.4 Closing the Unitarity Triangle 

Recently, several studies pointed out a possible tension in the SM between the value of sin 2/3 predicted 
from \ex\ and AM S / AMd, and its direct measurement from Sjj^k s [H> HI' HI' Hi- ^ n * ms sec ti° n ; 
we illustrate how CMM effects can remove this tension, and simultaneously account for a sizeable 
CP-violating phase in the B s system. 

Due to the particular sensitivity of \ck\ to new-physics effects, either 9 or 4>k must be very small. 
We will thus consider the two limits 9 = and 4>k = 0. For each case, we will compare the value of 
sin2/? extracted from Sj/^k s with its determination from \V US \, \V c b\, |e^|, and AM S / AMd, obtained 
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Figure 4: Constraints on 9 from B physics observables. Black (gray) points indicate allowed regions 
in Set 2 (Set 1) parameter space. The Erst four plots show individual three-sigma constraints from 
(a) AM d , (b) S J/lpKs , (c) AM d /AM s setting cp K = 0, (d) AM d /AM s setting cp K = vr/2. Plots (e) and 
(f) show the combined (a,b,c) and (a,b,d) constraints, respectively. In the case of Set 1, the three-sigma 
constraint from AM S has also been imposed, excluding points outside the 1.2 < |2^>.b s | < 2.4 range 
(recall that Set 2 is not affected by this constraint). Imposing further the constraint from <p s would 
remove the gray points with 2<fi b s < and the black points with —1.9 < 1§b s ^ —1.5 for sin# below 
0.15, see Fig. Qright. Finally, Set 2 (Set 1) points above the black (gray) horizontal line are excluded 
by AM K . 
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Figure 5: One-sigma constraints on the UT from Sjj^ Ks (light gray), \ex\ (gray), and AM S / AMd (dark 
gray) in the SM. The one-sigma region determined from \ V US \, \ V c b\, \zk\, and AM s /AMd assuming the 
SM is shown in black, and its shift due to CMM effects is indicated in dashed red. Left: Scenario I, 
(9 = 0, 4> Bs = 0.7. Right: Scenario II, 6 = 0.1, <j) Bs = 4>B d = 0.7. CMM inputs: Set 1. 



inverting the following expressions with respect to sin 2(5 and Rt- 



, , M K F%B K J G 2 F M$y 2 2 



V cb \ 2 R 2 t sm(2f3) m S (x t ) 



+ 2R t sm/5 (r)3So(x c ,x t ) - r)iS (x c )) 



2 sm 2fe sin (2(9) S™>(ri,r 3 ) >, 

8m| r 



(40) 



AM, ,Mg A/( A: i+ Xcos2 ^ s cos 2 0) 2 + (-2fe 2j R t sin ( 5|K s | 2 -^sin2^ s cos 2 r : 



2 



AM d 5 M Bd ^ (#2 cog 2{3\V US \ 2 + X cos 2cj) Bd sin 2 0) 2 + (i? 2 sin 2(5\V US \ 2 - X sin 2<j> Bd sin 2 1 

Here fci = 1 + |K S | 2 (1 - 2R t cos (5), k 2 = 1 + \V US \ 2 (1 - R t cos (5), 
7: 2 a 2 s (M s ) m S^( ri ,r z ) 

2\V cb \ 2 G 2 F M 2 ,m 2 -rS (x t y ^ 

and Rt = \VtdV^ b \/\V c dV* b \ is a side of the unitarity triangle (UT). The above expressions hold to 0.5% 
accuracy. In the SM, this leads to sin(2/? €A ') = 0.81 1qo9 with the inclusive \V c b\ determination of 
Tab. g, and to sin(2/? £A ') = 0.98 ±£11 if the exclusive determination from B -> D*tv decays, |V c fe| cxcl = 
(38.8 ± 1.1) • 10" 3 H|, is used instead. Note that \V cb \ incl does not lead to any significant deviation 



with respect to S^K Ks , while a tension is indeed observed with the smaller value |V^{,| excl . In order 
to illustrate how CMM effects can compensate for a low \Vd,\ input in UT analyses, we will adopt the 
averaged value of Ref. \ V cb\ LS = (41.0 ± 0.63) • 10~ 3 . In the following, we use the CMM input 
parameters of Set 1. All errors are treated as gaussian. 

6 = 0: CMM effects in R t Since for 9 = there are no effects in K and Bd mixing, CMM 
contributions enter the UT only via AM S . From Fig. |5|4eft, one sees that Rt has to increase in order to 
close the UT. This requires a CP-violating phase 2<^ Bs £ [1-2, 1.8], taking into account the three-sigma 
constraints on (pB s from AM S and <fi s . The dashed red curve shows Rt for 4>b s = 0.7, such that the UT 
determined from \ek\ and AM S / AMd agrees with the sin 2(5 measurement from Sj/^^s- 
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cf) K = 0, 6 = 0.1: CMM effects in R t and (3 In this second case, CMM effects affect both 
AM s /AM(j and Sj/^ Ks . For a fixed angle 9, the UT can be closed by adapting the CMM phase 
§B a — 4>B d - The resulting apex of the UT is shown by the intersection of the dashed red lines in 
Fig. ||-right for 9 = 0.1 and 4>b s = 0.7. For any value of 9 allowed by the constraints from B^ and B s 



observables in Sec. 5.3, one can find a phase (pB s to close the UT. 



Deviations from these two limit cases, i.e., small but nonzero 9 or (fix values, rapidly generate CMM 
effects in |e^| as well (Fig. |-left). These can lower the band from the \ex\ constraint in the (p, rj) plane, 
directly making up for the low |V^,| input value. 



6 Conclusions 

Grand-unified theories introduce relations among quark and lepton masses and mixings. Motivated by 
the large atmospheric mixing angle in the neutrino sector, several studies focussed on the consequences 
of the SU(5) Yukawa relation = Yj in b — > s transitions. In this work, we considered corrections to 
this relation which are essential to account for the observed light quark and lepton masses. In particular, 
we investigated the effects on s — ► d and b — ► d transitions of the additional rotation of the dn and sr 
quarks. This deviation with respect to the PMNS matrix, denoted by U, can be parameterized by an 
additional mixing angle 9 (see Eqs. (OJH))- 

In our analysis, we focussed on models with small Higgs representations; a modified version of the 
CMM model served as our specific scenario. In this setup, the differences between the down-quark 
and charged-lepton masses are naturally explained by dimension-five Yukawa operators. The associated 
supplementary rotation 9 was constrained from K — K and — B^ mixing observables. In particular, 
we found that, in the absence of fortuitous cancellations among the new phases in the matrix U, \ex\ 
sets a stringent bound on 9, # max ~ 0(1°). Consequently, in the basis where the charged-lepton Yukawa 
couplings are diagonal, the matrix DeY^ + Yj-D e + S^Yj-Yo- (in the notations of Eqs. (]5|,[lC|)) must be 
diagonal as well. Barring cancellations, this implies that the flavor structure of the couplings which 
modify the Yukawa unification must be similar to that of the initial terms. In other words, in the 
corrected relation Y^ = Yj + 5^ Y CT (Eq. (|Io|)), the three matrices Y CT , Y^, and Yj must be essentially 
aligned. Constraints from B physics observables (AM^, Sji^k 3 i and AM S / AM^) were also analyzed, 
and shown to imply the looser bound # max ~ 0(10°). 

While we have worked out this analysis for a specific GUT model, our results hold in general 
for models with small Higgs representations: large effects of the neutrino mixing angles on b^ — > sr 
transitions lead to large effects in 6r — > d^ and sr — > dn transitions for natural values of the parameters, 
once the mass relations for the light quarks and leptons are corrected. An efficient mechanism is naturally 
needed to render the mixing among right-handed <i-quarks visible. In the CMM model, this mechanism 
is provided by the fast SO(10) running of the dn soft mass matrix, which generates the large universality 
breaking Aj at the electroweak scale. Of course, other GUT scenarios could include additional sources 
of flavor and CP violation inducing effects in \ck\- These could soften the constraints on 9. Yet they 
would have to be fairly fine-tuned to cancel the potentially large impact of the corrections from the dn 
rotation matrix Rd (Eq. ([Tl|)). 

Interestingly, the correction operators which are of importance for proton decay but contribute 
equally to the fermion masses ought to have a different flavor structure in order to be in agreement 



with the experimental limit 17]. Both types of operators are generically present in GUTs. Hence, our 
analysis is an important step in establishing a consistent grand-unified model. 

Finally, we also considered the possible tension between the value of sin 2/3 predicted from \ex\ and 
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AM S / AM,; in the SM and its direct measurement from Sj/^k s > raised by the authors of Refs. [21, [?7|, 
|38| , 39]. We illustrated how CMM effects can remove this tension, and simultaneously reduce the 2.2a 
discrepancy observed recently in the B s — B s mixing phase. 
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Appendix: Loop Functions 



S (x c ) = x c , 

4x t — llx 2 + xl 3xf log(xt) 



So(x t ) 

S (x c ,x t ) = x 
F{x,y) = - 
G(x,y) 



4(1 -x t ) 2 2(1 -x t ) 3 ' 

3x t 3x 2 log x t 



log — 



x c 4(1 4(1 -x t f 



1 



1 



(x-l)(y-l) x-y 



xlnx 



y In y 



l 



+ 



1 



(x-l)(y-l) x-y 
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L o( x ,y) = j^G{x,y) - -F(x,y), 



(x - 1)2 (y - 1)2 
x 2 In x y 2 In y 



(x-l) 2 {y-lf 



5 (s) (a;,y) = L (x,x) - 2L (x,y) + L (y,y). 



(43) 
(44) 
(45) 
(46) 
(47) 
(48) 
(49) 
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